from pylab import *

seno=imread('seno_gray.png')[:,:,0]
coseno=imread('coseno_gray.png')[:,:,0]

senop=imread('seno_pared.png')[:,:,0]
cosenop=imread('coseno_pared.png')[:,:,0]

seno=2*(.5-seno)
coseno=2*(.5-coseno)


senop=2*(.5-senop)
cosenop=2*(.5-cosenop)



fase=arctan2(seno,coseno)
pase=arctan2(senop,cosenop)


figure(), imshow(fase),gray()
figure(), imshow(pase),gray()


dfase=fase-pase
dfase=where(dfase<0,dfase+2*pi,dfase)

figure(), imshow(dfase),gray()

